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Abstract: We compute the circular Wilson loop of A/" = 4 SYM theory at large N in 
the rank k symmetric and antisymmetric tensor representations. Using a quadratic 
Hermitian matrix model we obtain expressions for all values of the 't Hooft coupling. 
At large and small couplings we give explicit formulae and reproduce supergravity 
results from both D3 and D5 branes within a systematic framework. 



1. Background and motivation 

Wilson loops played a key role in the early development of the AdS/CFT correspon- 
dence [1]. The identification of the Wilson loop in A/" = 4 super Yang- Mills (SYM) 
theory with the dual fundamental string [2,3] connected naturally with ideas of un- 
derstanding the 'confining' strings of large N gauge theory as fundamental strings [4] . 
At finite temperature, thermal Wilson loops enabled the identification of a decon- 
finement transition in the theory [5] and allowed the computation of the thermal 
interquark potential and screening length [6-8]. 

There has recently been a renewed interest in Wilson loop operators in the 
AdS/CFT correspondence, with an underlying theme of understanding loops in 
higher representations of SU{N). This has seen a convergence of various differ- 
ent research directions. Firstly, circular Wilson loops are half BPS operators and 
have been conjectured to be computed exactly, up to instanton corrections [9], to all 
orders in N and A = Oym^ i from a quadratic Hermitian matrix model [10, 11]. This 
matrix model makes nontrivial predictions for Wilson loops in higher representations 
which constitute a test of the AdS/CFT correspondence incorporating higher genus 
string effects [12]. 

Recently, the matrix model computation of a A; winding circular Wilson loop 
was reproduced using a D3 brane embedded in AdS^ carrying k units of worldvol- 
ume fiux [13]. This beautiful match lead to further questions as it was not clear 
which rank k representation of SU{N) the D3 brane was computing. A complete 
dictionary was later proposed in [14] relating supersymmetric Wilson loops in an 
arbitrary representation to D3 and D5 branes in the dual geometry. According to 
this dictionary, a single D3 brane with k units of flux is dual to a Wilson loop in 
the rank k symmetric representation. It was therefore an open question to explain 
why the k winding computation of [13], a priori certainly not the same as the rank k 
symmetric respresentation, reproduced the supergravity result. This confusion lead 
us to the study of general symmetric representations in this paper. 

A second line of interest has originated from the description of certain half BPS 
operators of A/" = 4 SYM theory carrying angular momentum as 'bubbling geome- 
tries' [15]. Wilson lines provide another family of half BPS operators and so it is 
natural to search for a bubbling geometry formalism describing the backreaction of 
the D3 and D5 branes, as well as a field theory understanding in terms of fermion 
droplets. Exciting progress in this direction is currently being made [16-18]. This 
line of interest inspired the succesful matching of a bulk D5 brane computation with 
a matrix model computation of the rank k antisymmetric Wilson loop [19]. How- 
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ever, that matrix model computation required an assumption that is not a priori 
justified and does not appear to indicate corrections to the supergravity result. This 
provided the motivation for us to do an honest matrix model computation of the 
antisymmetric representation in this paper. 

Thirdly, computing the strong coupling thermal Wilson loops in higher repre- 
sentations is crucial for mapping out the phase structure of A/" = 4 SYM theory 
as a function of 't Hooft coupling and temperature [20]. It is an interesting open 
question whether there is a phase transition in the high temperature phase between 
weak and strong coupling physics, with implications for both black hole and plasma 
physics [21]. We recently initiated the required study of D3 and D5 branes in the 
AdS-Schwarzschild black hole background [22]. The results found for D5 branes in 
particular were curiously similar to those for the BPS circular spatial loop at zero 
temperature [19]. 

There have also been recent works on single Wilson loops in the fundamental 
representation [23, 24] . 

In this paper we study circular Wilson loops in rank k symmetric and anti- 
symmetric representations. Our computations will be in the large N limit of the 
Hermitian matrix model of [10,11]. The expressions we derive will be exact in A; 
we will take the large A limit only to obtain analytic results that may be compared 
with supergravity computations. A key formula in our framework is the contour 
integral representation of the symmetric and antisymmetric traces ( |2.26| ), which we 
then evaluate at large using a saddle point approximation. We go on to reproduce 
the known supergravity results for D3 and D5 branes [13, 19] as specific limits of this 
integral. 

The results in this paper provide a check of the proposed holographic dictionary 
for general representations [14], and furthermore allows one to compute the correc- 
tions to the supergravity limit. This last fact clarifies the regime of validity of the 
supergravity computations. In the case of symmetric representations, we describe an 
interesting phenomenon whereby the nature of the saddle point changes qualitatively 
as a function oik/N and A. We finish by deriving an explicit formula for the Wilson 
loops at weak coupling. 

Note: As we were typing up this work, the preprint [18] appeared which has 
some overlap with our results. 
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2. Circular Wilson loops 



The circular loop in the fundamental representation was first computed by Erickson, 
Semenoff and Zarembo [10] in perturbation theory by resumming a specific class 
of planar rainbow and ladder diagrams. It was further shown by these authors 
that cancellation of coordinate dependence in the planar Wick contractions makes it 
possible to map the result onto the solution of a large N Hermitian matrix model with 
quadratic potential. Subsequently this result was argued by Drukker and Gross [11] 
to be the exact answer, up to nonperturbative instanton contributions [9], to all 
orders in A = Oym-^ ^"^^ l/N. 

In this section we calculate BPS circular Wilson loops in the k-th rank antisym- 
metric and symmetric tensor representations of the SU(A^) A/" = 4 theory. Specifi- 
cally, we are interested in the expectation value of the operator 



[2.1] 



where TZ denotes the representation and 9^ is a constant unit vector in M^. This 
choice breaks the 5*0(6) R-symmetry of the J\f = 4 theory down to 5*0(5). The 
dimensions of the rank k symmetric and antisymmetric representions are respectively 
dim[5'fc] = ^l^j^j^I^y^ and dim[y4fe] = ^T(^ifc)T- The contour C will be chosen to be a 
spatial circle at zero temperature; this is then a BPS Wilson loop. 

At finite temperature, the above operator with the contour C chosen to be the 
Euclidean thermal circle yields the Maldacena-Polyakov loop. This is known to be an 
order parameter for the confinement /deconfinement transition in the strongly cou- 
pled field theory [5]. The Maldacena-Polyakov loop is not BPS in any natural sense 
since the finite temperature theory is not supersymmetric. Hence one does not expect 
the Maldacena-Polyakov loop to be related to the Hermitian matrix integral which 
computes zero temperature circular BPS Wilson loops. However, as mentioned in 
the introduction, the works of [22] and [19] suggest a tantalising connection. In par- 
ticular, the operator ( p.l| ) was computed in [22] using a probe D5 brane in the AdS- 



Schwarzschild black hole background yielding the Maldacena-Polyakov loop for the 
finite temperature theory in the antisymmetric tensor representation. Remarkably, 
up to an overall temperature dependent normalisation factor, the result was found 
to be exactly the same as the zero temperature circular BPS Wilson loop in the anti- 
symmetric tensor representation [19]. This suggests that at least in the supergravity 
(A = Qym-^ ^ 1) approximation there is a connection between Maldacena-Polyakov 
loops and Hermitian matrix models. This paper, however, will be concerned with 
the zero temperature theory and supersymmetric loops. 
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2.1 Circular Wilson loops from the matrix model 



The circular supersymmetric Wilson loop in the fundamental representation of SU(A^) 
was calculated in [10] and [11] and shown to be computed by a Gaussian Hermitian 
matrix model. The latter is defined via the partition function 



J[dM] 



exp ( -^Tr[M2] ) , (2.2) 



where M is a,n N x N Hermitian matrix and A = Qym^ is the 't Hooft couphng of 
the = 4 gauge theory. 

Assuming the conjectured relation between BPS Wilson loops and the matrix 
model for generic representations, circular Wilson loops in any representation TZ of 
the SU (N) gauge theory are then given by 

^ (Tr7e[t/]) = (Tr^^e^)^^. (2.3) 



dim[n] ' " dim[7^] 

The expectation value on the right hand side is computed in the Gaussian matrix 
model above. This is a remarkable conjecture which we will employ to extract the 
Wilson loops in symmetric and antisymmetric tensor representations. 

In an eigenvalue basis where M = diagjmi, m2, ms, . . . , mjv}, the expectation 
values of the operators in the completely antisymmetric tensor representation can be 
expressed as 

{Tr Ak[e^])mm = (expfmii + + . . . + mii^])mm = {0ik)mm , (2.4) 

l<n<i2<---<ife<A'' 

which are simply /c-th order symmetric polynomials a.^ in the eigenvalues {e"^'}. The 
{oik} are encapsulated in the natural generating function for symmetric polynomials 
which is the spectral curve or characteristic polynomial associated to the matrix 

N 

FA{t) = det{t + e^) = Hit + e'"^. (2.5) 

i=l 

It is easily seen that is the coefficient of t^^'^' in the spectral curve polynomial. It 
is most useful to express this relation as a contour integral in the complex t plane 

At any finite N the contour may be centred either around the origin or the point at 
infinity. Both definitions yield the same answer up to an overall sign. For infinite N 
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we will see that it is more convenient to think of the contour as encircling the point 
at infinity. 

The A;-th rank symmetric representations are given by the polynomials 0"^ in the 
eigenvalues {e"**} 

(TrSfc [e^])mm = X] (exp[mii + rui^ + . . . + J)^^ = {ak)mm ■ (2.7) 

i<n<«2<-<jfc<JV 

The generating functional for these polynomials is just the inverse of the character- 
istic polynomial^. We find it convenient to define it as 

1 ^ 1 

As in the antisymmetric case, we can extract the polynomials by a contour integral 
around the origin 

Note that the generating function Fs(t) decays as for large \t\, and has simple 
poles given by the positions of the eigenvalues of the matrix e^^ , namely at t = e~™*. 
Near the origin Fs{t) ~ 1 + ait + a^t"^ + . . ., and hence the Ok are directly obtained 
from the above contour integral. 

Our strategy will be to compute the expectation value of the generating functions 
in the matrix model 



{FA,s{t))mm — 




(2.10) 



where, as usual we have written the matrix integral in the eigenvalue basis at the 
expense of introducing a Jacobian factor, namely the Vandermonde determinant 

A^= n {rn,-m,f. (2.11) 

l<i<j<N 

In principle it should be possible to evaluate the above integral for any finite N and 
A by the method of orthogonal polynomials and extract the expectation values of 
BPS Wilson loops. However, in this paper we explore these expressions in the limit 
of large N and fixed A. which is relatively simpler and already very interesting. 

^We would like to thank Asad Naqvi for pointing this out to us. 
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2.2 The large limit 



To take the large N limit of the multi dimensional integral ( |2.10|) , we first rewrite it 
so as to make the dependence of various terms explicit 



1 r ^ 

^ / ll[dmi] expi-SA,si{m^})) , (2.12) 
^ i=i 



where 

Sa,s 



2N 

X 



N N ^ 

i=l i<j i=l ^ 



-log(t + e'"0 for Sa 
+ log(l -te™') for Ss 



(2.13) 



Clearly the tree level quadratic term is of order A^^, as is the Vandermonde repulsive 
interaction between the eigenvalues. On the other hand the insertion of the char- 
acteristic polynomial term in the partition function only contributes at order A^ in 
Sa,s- Therefore, in the large A^ limit, the integral will be dominated by the usual 
Gaussian model saddle point configuration wherein the eigenvalues condense onto 
a cut in the complex eigenvalue plane with a spectral density characterised by the 
Wigner semi-circle distribution (see e.g. [25] and references therein) 



p{m) 



7tX 



(2.14) 



The spectral density is normalised so that 



dm p{m) 



1. 



(2.15) 



The large A^ expectation values of the generating functionals are thus 
{FAit)Um = (exp(Trlog[t + e^']))^^ 



exp 



A^ 



dx p{x) \og{t + e^) 



(2.16) 



and 



(exp(-Trlog[l-te^^]))^^ 

-A^ / dx p{x) log(l - te" 



exp 



(2.17) 



At this point we note that in the strict large A^ limit there is a drastic change 
in the analytic structure of the characteristic polynomial {FA{t))mm- This can be 
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understood in two different ways. Firstly, for any finite N the representation of 
tlie polynomial as exp(Trlog(t + e^)) yields only an A^-th order polynomial upon a 
formal power series expansion around the origin. The coefficients of higher powers of 
if: cancel due to trace relations for large powers of the N x N matrix e^. At infinite 
A^, as is evident from ( p.l6| ), the power series in t around the origin has infinitely 
many terms. This is due to the fact that the N ^ oo limit has already been taken 
and one does not see traces involving powers of e*^ that are larger than N. 

In addition, we note that on the real axis for t < — the expression develops 
an imaginary part. This can be traced to a finite branch cut discontinuity along 
the negative real axis. The discrete zeros of the characteristic polynomial, which are 
the eigenvalues of , coalesce into a continuum in the large N limit producing the 
branch cut along — ^ t < —e^^. Hence at infinite our expression for the 
spectral curve does not yield a simple polynomial since it has a branch cut singularity 
in the t plane along with a pole of order N at infinity. 

The generating function {Fs(t))mm for the symmetric Wilson loops also has a 
branch cut on the real axis for < ^ < e^, which is due to the coalescing of 
simple poles in the large A^ limit. 

The large A^ circular Wilson loop in the rank k antisymmetric representation is 
then 



(TrAj^]) ^ k\{N-k)\ r dt 
dim[Ak] N\ Jc 



i — t^-^ exp A^ / dx p{x)\og{l + t-^' 
Jc 2vrz 



(2.18) 



where we have absorbed the factor of t^, see ( |2.6| ), into the exponent and the contour 
of integration C is to be specified below. For the rank k symmetric loops at large 
A^, we have similarly 



(Tr^Jf/]) k\{N-l)\ / dt 1 



dim[5fc] (A^ + A; - 1)! Jc 2m t^+i 



-A^ / p(x) log(l - t) 



(2.19) 



The problem of computing the Wilson loops is reduced to computing residues. 
This is quite straightforward in principle for finite k = 1,2,... A more interesting 
question is whether it provides a useful way to explore the behaviour of these loops 
for generic values of k that are comparable to A^. 

2.3 Wilson loops at A; ~ A^ 

A compact expression for the rank k Wilson loops is obtained after a t 1/t variable 
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change in ( p.l8 



) so that for the symmetric (S) and antisymmetric (A) loops we have 



{TTs„AM) = ^^£dt-^ exp tN j ^dxp{x)\og{lTe''t) 



(2.20) 



where the — sign corresponds to the symmetric representation. The contour C is 
chosen around t = and this encapsulates the prescription that the /c- Wilson loops 
are the coefficients of positive powers of t in the Laurent series around t = 0. Indeed, 
this prescription works as we obtain 



{TtsuaAU]) = N / dxp{x)e^ = {Tt[U]) 



{Tts,,aAU]) = ^(Trf/)2 ± l-N dxp{x)e^^ 



(2.21) 



(2.22) 



- {{TiUy ± (Trf/2)) . (2.23) 



This is not surprising but provides a simple check of the prescription. Note that 
strictly speaking, the expressions ( |2.21| ) to ( |2.23| ) are not quite reliable within the 
large limit. We have kept terms which scale differently with N, while replacing 
((Trt/)") with (Trf/)" which is of course a consequence of the N —>■ oo limit. In 
any case, the direct approach of evaluating explicitly the residue at the origin is not 
useful for computing Wilson loops with k comparable to A^. To access these we need 
to evaluate the saddle point of the integrand at large A^. 



2.4 Large A^ saddle point 



The integral ( |2.2U| ) yields the Wilson loops in the large A^ theory and therefore only 



strictly makes sense in that limit. In the N ^ oo limit we only need to consider the 
contribution from the saddle points of the integrand. 

We are interested in the regime — »• oo with k/N fixed. We introduce the 
fraction 

(2.24) 

and treat it as a continuous variable in the large A^ limit. In addition, it is natural 
to perform the variable change 

t = e^\ (2.25) 

which takes us from the complex plane to the cylinder. The points t = and t = oo 
are mapped to the points z = ^oo respectively. 
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For the antisymmetric case the branch cut singularity lying along the interval 
(— e^, — e~^) in the t plane is mapped to the line joining z = —1 + i-^ and 
z = 1 + and images thereof under translations by multiples of in the z plane. 

In the symmetric case the branch cut along < t < is mapped to the interval 
— 1 < 2; < +1 along the real axis and its images in the z plane. 

Using the explicit form of the Wigner semi-circle distribution and appropriately 
rescaling variables we find from ( p. 20 ) 



[Tts^aAU]) = ^,j^dz exp tN Q j'^ dx v/r^log(l T e-^(^-^)) ± fVXz^ 

(2.26) 

The contour C lies to the left of the branch cut, in the region Re 2; < — 1 and winds 
once around the z plane cylinder. 

In the N —>■ 00 limit, we expect the integral to be dominated by the saddle points 
of the exponent. Hence we will first look for these. The saddle point equations are 

- dx ~ , ± / = . (2.27) 

The choice of upper signs yields the symmetric case. These equations will in general 
have solutions in the complex z plane. Note that for certain values of z we need to 
worry about potential divergences in the integrand in (|2.27| ). The integral is then 
given by its Cauchy principal value along with an imaginary part. Such situations 
occur precisely when the solution lies on the branch cut and is not a stationary 
phase point. An explicit analysis of these equations is hard for generic A. Instead 
we first study the limit of large A which can be directly compared with supergravity 
calculations. 



2.5 Supergravity limit 



The supergravity limit, A ^ 00, leads to a drastic simplification of the exponent in 
( |2.26|) and the resulting saddle point equations (|2.27| ). For the rest of this section. 



we will treat the antisymmetric and symmetric cases separately. 



Antisymmetric: This corresponds to the equation ( p.27| ) with the lower sign 



choice. Separating out the saddle point equation into real and imaginary pieces, we 
find that the imaginary part of the equation leads to the condition 



sm 



v^Im;z) = 0. (2.28) 
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While this has the general solution Im z = nir, we will see below that only the 
solution Im z = (and images) is picked out as the saddle point. 



Now we turn to the real part of the saddle point equation. We first divide the 
z plane into three regions: Re2<— 1,— l<Re2<l and Re ^ > 1. In each of 
these three regions the limit of large A leads to distinct asymptotic expressions. The 
saddle point equation can be expanded in powers of in the large A regime to 
yield 

-/ + 0(e-^(l^''"l-^)) = [Rez<-1], (2.29) 
- dxyr^^-/ +0(e~^(l^""l-^)) = [-l<Rez<l], (2.30) 

J-1 

and finally 

(1 - /) + 0(e-^(l^^ ^1-^)) =0 [Re ^ > 1] . (2.31) 

It is clear that in the regions Re 2; > 1 and Re 2; < 1 there exist no consistent 
solutions in the strict A ^ 00 limit for fixed /. 



On the other hand, in the region — 1 < Re z < 1 we have a nontrivial equation 
that in general has solutions. Recalling that the branch cut extends between z = 
1 + i-^, we conclude that the stationary phase point hes at 



— 1 + and z 



Im z = 0, on the real axis at 



z e 



where z satisfies 



cos ^ z 



(2.32) 
(2.33) 



SVl -P= 7r(l-/), 
It is natural to rewrite this in an angular parametrisation 

z = cosek, (2.34) 

where the subscript k associates the solution to the rank k representation. This gives 

7i{l~ f) = i9k-sm9kCos9k). (2.35) 

Note that under k ^ N — k, or equivalently / — 1 — /, the angular variable 6k trans- 
forms to IT — 6k- We will see in a moment that under this transformation the integral 
( p.26| ) remains invariant, reflecting the fact that antisymmetric representations of 
rank k and N — k are in fact identical. 



The saddle point value of the integral is determined by the exponent in ( ]2.2(j| ). 
Putting z = z and taking the large A limit of ( p. 26 ) we obtain that the antisymmetric 
Wilson loop at large k N evaluates to 



dim[Ak 



■(Tr^Jf/])=exp 



2N 
37 



VA sin^ 6k 



(2.36) 
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Here we have only kept the leading terms in the strict large A limit in the exponent. 
There are several noteworthy points about this result. The first is that it matches 
the answer obtained by a D5 probe brane calculation in AdS^ x 5*^ [19] and more 
surprisingly, the D5 brane probe calculation of the Polyakov loop in the AdS black 
hole background [22]. Although this answer was also obtained by a matrix model 
calculation in [19], the logic of that calculation is unclear ^. In particular, our result 
will have corrections involving powers of in the exponent. Such corrections 

were also present in the matrix model computation of [11] for the Wilson loop in the 
fundamental representation. Note also that the result is invariant under 6k —>■ tc — 6k 
reflecting the k —>■ N — k symmetry for completely antisymmetric representations. 
Finally, as noted in [22] and [19], the above form of the action implies terms that 
appear nonanalytic from the point of view of the 1/N expansion. Specifically, one 
obtains contributions involving powers of {j^)'^^^ upon expanding the saddle point 
exponent in powers oi f = k/N . 

Symmetric: For the rank k symmetric representation, the saddle point equation 
corresponds to choosing the upper signs in ( |2.27| ). This choice of sign completely 
changes the nature of solutions. First, by separating the equation into its real and 
imaginary parts we may show that there are no saddle points away from the real 
axis on the z plane. Recall that in the symmetric case the branch cut is along the 
interval —l<z<l. 



If we start with A ^ 1, then we find a saddle point on the real z axis with 
Re 2; < — 1. This saddle point can be treated analytically and we will do so in the 
following section. As we increase A, however, we find that the saddle point moves 
in the positive direction along the real axis. At a critical value, A = Ac, the saddle 
point hits the branch point at Re z = — 1 and appears to move off into the second 
sheet. For a given /, from ( p.27| ) we have that the critical point is reached when 



(2.37) 

Solving this equation we obtain a curve that splits the (/, A) plane into regions where 
a saddle point exists for negative real z and where it does not. This curve is shown 
in figure 1. 



^The calculation presented in [19] was in fact a computation of the cumulative distribution of 
the first k eigenvalues in the cut using the large N distribution ^ ^'^^ dxx\J\ — . Although 
this matches the correct large A result, it differs from the correct finite A result as we already see 
from the 0(e~^) corrections in our approach. Obtaining the correct answer using the approach 
of [19] involves computing the cumulative expectation value of k distinct eigenvalues sprinkled ran- 
domly along the distribution governed by the Hermitian matrix model. These will give appropriate 
corrections to the result of [19] 
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Figure 1: The existence of a saddle point at negative real z depends on / and A as 
shown. 

The asymptotics of the curve in figure 1 may be derived analytically. We have 

2 

/ - + • ■ • as Ac ^ , 

Ac 





No Saddle 


Saddle 





In the region where the saddle point exists, we can compute the location of 



the saddle point numerically and evaluate (|2.26| ) at the saddle, thereby giving us 
the Wilson loop expectation value. However, at any fixed /, once the coupling is 
sufficiently strong there are no saddle points on the first sheet of the complex plane. 
We are aiming for analytic expressions in the \ ^ oo limit. Since it is not clear how 
to pick up saddle point contributions from the second Riemann sheet, we will adopt 
a different methodology to evaluate ( p.26| ) for the symmetric representation. 



The integral has no singularitites at Re z +oo (corresponding to t ^ cxd in 
our initial coordinates). Its only singularity is the branch cut on the real axis. We 
may therefore deform the contour to wrap the branch cut so that the integral reduces 
to the discontinuity across the cut 



(Trsjf/]) = ^Im y dy exp -A^ dx Vl - log(e^^ 



) + 



9 r 

+ - J dx v/r^^log(e^^' - e^y) + 2i J 



dx\f\ 



X 



2 + f^fX^ 



(2.39) 



The Wilson loop is thus given by a definite integral over the real variable y and we 
can unambiguously evaluate it using the stationary phase method at large A^, by 
looking for the saddle points of the integrand analytically continued to the complex 
plane. 
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We find two sets of saddle points in tlie large A limit, one whose location depends 
only on the combination f\/X in the strong coupling limit, and another set whose 
location is independent of A in this limit. 



The first of these saddle points occurs along the negative real axis and solves the 
equation 

7T J _i ^y^y — gvAx 

The factor of 4 in this expression arises for the following reason. Analytically con- 
tinuing the exponent of ( |2.39| ) into the complex y plane and in particular to large 



negative values of y yields an extra imaginary contribution from the logarithm which 
then doubles the coefficient of the third integral in the exponent. 

At large A and a fixed value of where 

(2.41) 



we find a solution to the saddle point equation ( |2.40| ) along the negative real axis at 

y = yi = -v/rr^+0(e-^(l^^l-i)). (2.42) 



Plugging this back into ( |2.39|) we obtain the contribution to the Wilson loop from 
this saddle point as 



Wi = exp 



2N (^kVI + ^2 + sinh"^ , (2.43) 



which is exactly the result found by Drukker and Fiol [13] using a probe D3 brane 
in ^^5*5 X S^. We have only kept the leading terms in the large A^, large A limit 
in the exponent. In particular we have ignored phases. This is only one of the two 
contributions to the symmetric Wilson loop. 

The integrand ( p.39| ) appears to have yet another saddle point at large A in 
the complex plane, whose position is independent of A in this limit. To see this we 
directly take the large A limit of the exponent (|2.39| ). Once we have taken this limit 



we can look for the saddle point. The condition turns out to be very similar to ( |2.30| ), 
obtained earlier in the antisymmetric case 



TT 



dxVl - + / = , (2.44) 

-1 



except that y can be anywhere in the complex plane. It is obvious that there are 
no real solutions to this equation. This is a key difference between the symmetric 
and antisymmetric representations. The expressions obtained are only superficially 
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similar. It may be checked numerically that the equation has solutions in the complex 
plane. In terms of a comj)/ea: angular variable, the complex saddle point in the y plane 



y2 = cos0fcGC, (2.45) 

is determined by 

(1 + /)7r = 0fc - sin0fccos0fc- (2.46) 

This does not have any solutions for real ipk between and vr. The main difference 
from the antisymmetric case ( ^.351 ) is that now the left hand side is greater than 
unity. 



This solution contributes to the symmetric Wilson loop in the limit of large A as 

2N 



W2 = exp 



VX Re(sin^ 0^) 



(2.47) 



where we have ignored phases and terms that are subleading in the exponent in the 
large A regime. We have checked numerically that for generic values of A and / the 
exponent turns out to be positive at the saddle point and hence is nonzero in the 
large N limit. 



The Wilson loop in the symmetric representation is given by the sum of these 
two stationary point values 

= dh^ ^^"'^ f^^^ = 1^1 + . (2.48) 

Depending on the specific values of f = k/N and A, one of the two stationary points 
provides the answer in the large N theory. The limit that is relevant for direct 
comparison with the D3 brane computation in supergravity is A — >■ 00 with k fixed. 
In this limit the saddle point Wi is dominant and therefore our computations match 
the D3 brane result of [13]. Another limit in which Wi dominates is the large A limit 
with k/N fixed, since its exponent scales as 2Nk'^ ~ A^A while the exponent of W2 
only scales as N\/X. For general values of A and /, either saddle could dominate. 



2.6 Perturbative limit 



Interesting explicit results may also be obtained in the opposite limit, A 0. In this 
regime, a unified treatment of the symmetric and antisymmetric cases is possible. 



The first step is to solve the saddle point condition ( |2.27| ) perturbatively in A 



As in the previous subsection, the saddle point must lie on the real axis: z = z & 
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One finds the following expansion to the lowest two orders 



1 



1 + 2/- 
/ 



+ 



(2.49) 



As elsewhere, the upper choice of signs corresponds to the symmetric case. 



Substituting the expansion ( p.49| ) into ( p.26| ) and expanding the expression in the 
exponent, we obtain the following result for the perturbative Wilson loop to leading 
order as A — 



W: 



[U] = exp 



(2.50) 



We can recognise the exponent of this result as the quadratic Casimir of the symmet- 
ric and antisymmetric representations of SU{N) in the large N limit with f = k/N 
kept fixed. This observation strongly suggests the following perturbative result for a 
general representation 



Wn[U]=exp 



(2.51) 



The appearance of the quadratic Casimir is not surprising and the linear dependence 
on A is consistent with a one-loop perturbative gauge theory result. The coefficient 
is a nontrivial result of the A/" = 4 gauge theory. 



A more detailed study at weak coupling in the symmetric case provides a pleasing 
consistency check of the strong coupling results. There exists a saddle point which 
is the weak coupling counterpart of the saddle ( p.42| ) that we matched with the 
supergravity D3 brane result. This may be seen most easily in the limit A with 
fixed K 1. In this case the saddle point equation (|2.40| ) is solved by 



y = y^ = -ViT^+0{k-^] 



(2.52) 



Evaluating the Wilson loop on this saddle point yields precisely the large k limit of 
the result ( p.43|) . Thus although we might appear to be far from the supergravity 
regime, we can still make contact with the D3 brane result of [13]. 



3. Conclusions 



In this paper we have used the quadratic Hermitian matrix model of [10,11] to com- 
pute circular higher rank Wilson loops in SU{N) Af = 4 SYM theory. In particular, 
we gave an expression for the rank k symmetric and antisymmetric representations 
in terms of a single contour integral. This integral expression is valid for all val- 
ues of the 't Hooft coupling and at large A^ may be computed in the saddle point 
approximation. 
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Taking the strong coupling limit, A — cxo, we reproduced analytically both 
the supergravity D3 brane result of Drukker and Fiol [13], using the symmetric 
representations, and the supergravity D5 brane result of Yamaguchi [19], using the 
antisymmetric representations. This provides nontrivial evidence for the Wilson loop 
dictionary proposed in [14] and of the AdS/CFT correspondence more generally: the 
match involves higher genus terms, as k/N was kept fixed in the large N limit. Our 
results also clear up confusions in the literature, in which the D3 brane result was 
initially matched with a multiply wound loop while the D5 brane result was matched 
with a matrix model computation whose validity was not a priori obvious. 

Away from the strong coupling limit we find corrections to the supergravity 
results governed by e^^. These corrections are computable, if unilluminating at 
present, within our framework. We also computed the Wilson loops in the small A 
limit, and were led to a simple formula that we expect to be true for an arbitrary 
tensor representation with mixed symmetries. 

The interpolation from weak to strong coupling appears to be uneventful in 
the antisymmetric case. In contrast, there is a very interesting phenomenon in the 
symmetric case. At a critical value of the coupling, the saddle point of the contour 
integral hits a branch point and moves off the real axis and into a second Riemann 
sheet. We have shown the dependence of the critical coupling on k/N, but it seems 
likely that there is interesting physics around this point still to be uncovered. 

The formalism we have presented can also be used at finite A^, with the ma- 
trix model expectation values evaluated using orthogonal polynomials rather than 
contour integrals. This may be another interesting direction for future work. 
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